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Summary 

In  the  first  four  months  of  the  Project,  various  aspects  of  light  self-localization  in 
reorientational  Kerr-like  media  such  as  nematic  liquid  crystals  have  been  examined,  including 
issues  related  to  the  propagation  of  narrow  beams  besides  the  paraxial  approximation  as  well 
as  self-acceleration  of  self-confined  light  beams  in  Kerr  anisotropic  dielectrics. 

The  samples  were  provided  through  collaborations  the  University  of  Southampton  (UK),  the 
Institut  Non  Lineare  de  Nice  (Sophia  Antipolis,  France))  and  BEAM-Engineering  (Orlando, 
USA).  After  a  brief  introduction,  we  recall  the  basic  concepts  on  self-localized  propagation  of 
light  beams  in  media  with  a  reorientational  response;  then  we  illustrate  numerical  and 
theoretical  analyses  and  some  experimental  results  on  self-steered  self-confined  light  beams, 
i.e.  self-accelerating  spatial  solitons,  in  nematic  liquid  crystals  (LC). 

1  Introduction 

Coupling  and  routing  light  for  optical  telecommunications  and  on-chip  signal  processing 
is  currently  mostly  based  on  static  couplers  or  opto-electric  and  electro-optical  converters. 
Experiments  in  LC  have  proven  that  two-dimensional  solitons  are  "robust"  entities,  can 
survive  index  non-homogeneities  (lens-like,  interfaces  or  other  solitons),  can  confine  weaker 
signals  with  a  smooth  transverse  profile,  can  be  steered  and  routed  in  space  by  externally 
controlled  parameters,  including  temperature,  voltage,  extra  beams.1'3  Owing  to  self-healing 
upon  propagation,4  liquid  crystals  soliton-induced  waveguides  can  tolerate  alignment  and 
angular  offsets,  as  well  as  polarization  noise.  While  optical  spatial  solitons  reported  so  far  in 


liquid  crystals  have  been  shown  to  react  to  external  stimuli,  their  performance  is  far  from 
optimal:  their  response  time  can  be  larger  than  100ms  and  they  exhibit  transverse  fluctuations 
due  to  non-locality  and  the  liquid  phase  of  the  supporting  medium.  Overcoming  the  current 
limits  is  necessary  to  realize  competitive  and  all-optical  signal  interconnects.  Based  on  recent 
progress  spatial  optical  solitons  in  nematic  liquid  crystals,  spatially  stable  waveguides  induced 
by  self-trapped  light  beams  can  be  envisioned  as  a  novel  photonic  platform  for  dynamically 
re-addressable  broadband  optical  signal  interconnects. 

As  anisotropic  molecular  systems  are  amenable  to  the  definition  and  tuning  of  arbitrary 
refractive  index  distributions  by  external  voltages  and  light-induced  reorientation,  in  this 
Project  we  aim  at  verifying  models  for  the  exploitation  of  electro-optic,  as  well  as  opto-optic, 
control  of  soliton  waveguides  and  circuits  in  nematic  liquid  crystals,  including  light-valves 
with  photosensitive  layers.5  The  merging  of  non-local/electro-optic/non-linear  material 
properties  and  engineered/tunable  refractive  index  distributions  is  meant  to  naturally  eliminate 
some  drawbacks  of  ffee-space  propagation  and  offer  more  versatility  than  built-in  integrated 
optical  interconnects.  Moreover,  the  intrinsic  self-collimation  (guidance)  characteristic  of 
bright  solitons  will  reduce  signal  attenuation,  ensure  a  smooth  transverse  profile  ideal  for  in¬ 
fibre  coupling,  reduce  polarization  noise  and  allow  large  deflections  in  engineered  structures. 

Such  new  photonic  platform  and  interconnect  approach  requires  an  in-depth  analysis  of 
the  intricate  interplay  between  non-linear  light  propagation  in  non-local  media,  the  interaction 
of  light  with  boundaries,  refractive  index  profiles,  external  voltages  and  other  light  beams. 

2  Methods,  Assumptions,  and  Procedures 

The  Project  aims  at  demonstrating  the  feasibility  of  optical  and  all-optical  signal 
processing  for  network  reconfiguration  in  compact  microsystems  based  on  light  self¬ 
localization  and  power-controlled  self-steering  using  liquid  crystal  technology  and  liquid 
crystalline  materials,  specifically  nematic  liquid  crystals. 

The  two  main  routes  of  investigation  undertaken  during  the  first  four  months  were: 

-  1:  Light  propagation  and  localization  in  non-local  molecular  media:  we  focus  on  theory, 
simulation,  design,  realization  and  characterization  of  individual  waveguides  based  on  spatial 
optical  solitons  in  nematic  liquid  crystals,  i.e.  Nematicons.6  Beam  propagation  models  and 
software  will  be  used  to  analyze  and  design  the  required  refractive  index  distributions, 
modelling  the  effects  of  non-locality  and  the  role  of  tight  (self-)  confinement.  We  will  identify 
and  prepare  the  samples  to  be  assembled  in  two  main  geometries  -  namely  planar  cells  and 
planar  photoconductive  light  valves  -  for  optoelectronic  characterization.  The  basic  elements 
will  be  characterized  in  terms  of  excitation,  profile,  spatial  response  etc. 

-  2:  Steering  and  routing  with  optical  spatial  solitons:  we  focus  on  theory,  simulation, 
design,  realization  and  characterization  of  different  ways  for  transverse  self-acceleration  of 
Nematicons,  obtained  to  date  by: 

-  mutual  interaction  of  solitons  in  nonlocal  dielectrics  710 

-  voltage  control  u'14; 

-  light  control 15~18. 


3.  Results  and  Discussion 

3.1  Non  paraxial  beam  propagation 

Since  the  excitation  and  propagation  of  optical  spatial  solitons  in  nematic  liquid  crystals 
involves  beam  self-confinement  via  a  reorientational  nonlinearity  with  a  strong  nonlocal 


character,  nematicons  are  expected  to  breath  in  waist  as  they  propagate  and  self-fosus  to 
transverse  size  comparable  or  even  narrower  than  the  wavelength  of  light.  The  latter  feature 
demands  a  detailed  investigation  of  narrow  beam  propagation,  as  their  reduced  waist  fails  to 
comply  with  the  usually  adopted  scalar  or  paraxial  approximation  and  the  slowly  varying 
envelope  approximation.20'21 

In  the  study  of  paraxial  beams  in  isotropic  homogeneous  media,  scalar  Hermite-Gaussian 
beams  are  essential,  as  they  form  a  complete  set:  when  dealing  with  one-directional  (1+1)D 
propagation,  every  initial  profde  can  be  expressed  as  a  superposition  of  HG  eigen-beams. 
Focusing  on  a  (1+1)D  geometry,  we  demonstrate  that  the  fully  vectorial  problem  of  light 
beam  evolution  in  uniaxials,  when  excited  by  extraordinary  waves  undergoing  a  finite  walk- 
off,  can  be  exactly  reduced  to  a  single  scalar  Helmholtz  equation  for  the  magnetic  field,  even 
in  the  non-paraxial  limit. 

We  consider  a  homogeneous  non-magnetic  crystalline  dielectric  encompassing  a  single  optic 
axis  and  in  the  presence  of  a  harmonic  electromagnetic  perturbation,  with  constitutive 

relations  B  =  pH  and  D  =  8o£*  E.  We  name  s±  and  s //  the  eigenvalues  of  e  with  multiplicity 

two  and  one,  respectively,  and  define  the  birefringence  sa  =  £//  ~e±.  For  simplicity,  we  take 

the  optic  axis  lying  in  the  plane  yz  and  a  (1+1)D  geometr,  such  that,  two  independent 
solutions  exist  for  Maxwell's  equations:  ordinary  (-o))  and  extraordinary  (-e)  waves, 
respectively.  Light  propagation  is  governed  by  two  Helmholtz's  equations.  In  a  rotated 
reference  system  introduced  to  handle  e-waves,  o-solutions  can  be  found  by  duality,  i.e. 
exchanging  E  with  H  and  viceversa  from  the  e-case  and  syz=0  (no  walkoff). 

Adopting  a  reference  system  xts,  with  s  along  the  Poynting  vector  S  and  t  orthogonal  to  it,  the 
electric  field  components  take  a  relatively  simple  form  and  the  flux  of  S  through  a  close 
volume  is  purely  reactive,  in  agreement  with  energy  conservation  in  the  stationary  regime.  For 
a  finite  beam  and  an  arbitrary  direction  z,  we  can  address  the  effects  of  anisotropy  on  the 
symmetry  of  the  Poynting  vector.  We  consider  an  input  magnetic  field  A  of  given  parity:  such 
solution  of  a  Helmholtz's  equation  will  retain  its  symmetry  while  propagating.  For  a  light 
beam  with  mean  wavevector  along  z  and  its  corresponding  Poynting  vector  along  s,  a 
symmetry  breaking  takes  place  in  anisotropic  dielectrics  where  the  walkoff  is  nonzero  and  its 
effects  increase  with  the  beam  divergenc,  i.e.  its  spatial  spectrum. 

Considering  a  narrow  beam  with  an  initial  profile  given  by  a  HG  function,  we  use  a  numerical 
expansion  of  the  angular  spectrum  to  solve  the  propagation  problem  exactly,  focusing  on  the 
lowest-order  Gaussian  input  profile  with  initial  waist  wo  and  a  flat  phase  profile. 
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Fig.  1  (a,  d)  Beam  waist  vs  z;  blue  solid  lines  and  red  symbols  refer  to  paraxial  Gaussian  beams  and  numerical 
solutions,  respectively,  (b,  e)  Hx2  vs  y/X,  normalized  to  its  peak  and  computed  for  (from  left  to  right)  z/  X  =0,  53, 
105,  158;  solid  and  dashed  lines  are  FIG  profile  and  numerical  solutions,  respectively,  (c,  f)  Relative  error 
between  the  profiles  in  (b,  e).  Top  and  bottom  graphs  are  computed  for  w0/X  =0.5  and  w0/  X  =2,  respectively. 
Flere  en  =2.89,  s±=2.25.  Inset:  unit  vectors  s  and  t. 


The  numerical  solutions  confirm  the  theoretical  predictions:  the  magnetic  field  walks-off  and 
retains  its  even  symmetry  even  in  the  highly  non-paraxial  limit.  As  visible  in  Fig.l,  going 
from  large  waists  down,  the  effects  of  the  second  partial  derivative  of  A  along  the  propagation 
z  are  negligible  up  to  the  wavelength  X;  for  w/X  =2  the  relative  error  is  <1\%,  whereas  for  w/X 
=0.5  the  differences  between  the  two  models  remain  small,  with  a  relative  error  on  waist 
<8\%  (Fig.l).  The  trend  is  confirmed  by  Fig.l(c-f)  plotting  the  relative  error  between  paraxial 
and  non-paraxial  Hx.  The  agreement  between  numerical  simulations  and  HG  profile  for  large 
w/X  confirms  that  the  slowly-varying-enevelope  approximation  needs  be  applied  along  in  the 
moving  reference  system.  Fig.  2(a)  plots  Re(S)  along  s,  corresponding  to  Hx  in  Fig.  1(a):  the 
energy  propagates  at  the  walk-off  angle.  Both  the  ratios  between  the  longitudinal  Es  and  the 
transverse  Et  electric  field  components  and  between  odd  and  even  parts  of  Re(S)  are  inversely 
proportional  to  wo  [Fig.2  (c,  d)],  with  the  former  one  order  of  magnitude  larger  than  the  latter. 
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Fig.  2  (a)  Computed  transverse  profile  of  Re(S)  along  s  and  (b)  relative  error  with  respect  to  the  paraxial 
expression  for  w A,  =0.5;  Re(S)  is  normalized  to  its  peak  for  fixed  z  and  each  curve  corresponds  to  the  same  z /X 
of  Fig.  1.  (c)  Ratio  of  maxima  of  |Et|  and  |ES  |  for  fixed  z  and  (d)  ratio  between  maxima  of  odd  and  even  parts  of 
Re(S)  vs  Wo/A. 


A  direct  comparison  of  Fig.  1  (c,  f)  with  Fig.  2  (c)  reveals  that  the  differences  between  scalar 
(paraxial)  and  full  vectorial  model  solutions  are  larger  than  the  corrections  due  to  the  second 
longitudinal  derivative,  in  agreement  with  Ref.  20. 

3.2  Spatial  optical  solitons  in  nematic  liquid  crystals 

In  many  cases  the  nonlinear  effects  are  related  to  the  anisotropy  of  the  material,  such  as  in 
reorientational  media  and  nematic  liquid  crystals.  ’  ’  ’  Spatial  solitons  in  the  paraxial 
regime  are  often  modelled  by  a  scalar  NL  Schroedinger  equation. 

We  extended  the  anlysis  summarized  in  the  previous  section  and  investigated  bright  single¬ 
hump  symmetrically-shaped  spatial  solitons  in  uniaxial  media  by  solving  Maxwell's 
equations  in  geometries  where  ordinary  and  extraordinary  polarization  components  are 
decoupled,  studying  non-perturbative  regimes  where  linear  and  NL  polarizations  become 
comparable  in  magnitude  and  showing  how  non-paraxiality  results  in  an  effective  nonlocality 
even  in  local  media. 

In  particular,  in  the  highly  NL  regime  relevant  to  nematic  liquid  crystals  (NLC)  -among  other 
NL  dielectrics-,  we  expect  soliton  self-steering  via  power-controlled  changes  in  walk-off. 

NLC  are  known  to  exhibit  very  high  effective  nonlinearities  (typically  106  cm2/W)  due  to  the 
particular  nature  of  light-matter  interactions.22  As  a  consequence,  in  uniaxial  NLC  the  optic 
axis  (or  molecular  director)  can  reorientate  in  space  by  reacting  to  the  electric  fields  of  light 
beams  in  the  milliwatt  or  even  microwatt  range17:  when  the  induced  rotation  becomes 
comparable  with  the  initial  distribution,  the  nonlinear  polarization  is  no  longer  negligible  with 
respect  to  the  linear  one.  Light  self-steering  owing  to  nonlinear  changes  in  walk-off  was 
experimentally  demonstrated  in  the  modulational  instability  regime.  23,  24  NLC  behave  as  a 


saturable  nonlocal  nonlinear  medium  able  to  support  stable  (2+l)D  spatial  solitons  at 
relatively  low  cw  excitations.  At  low  powers,  the  soliton  trajectories  do  not  bend  because  of 
the  negligible  walk-off  variations  associated  with  the  optically- induced  reorientation;11  at 
higher  powers,  beam  self-bending  can  be  expected. 

We  refer  to  a  planar  NLC  cell,  infinitely  extended  along  y  and  z  and  with  a  finite  thickness  L 
across  x.8  The  boundary  conditions  induce  a  homogeneous  director  distribution  in  the  absence 
of  excitation,  with  the  optic  axis  lying  in  the  plane  yz  at  an  angle  0o  with  respect  to  z. 
Maxwell’s  equations  need  be  solved  in  conjunction  with  those  ruling  the  director  distribution; 
since  the  electric  field  vector  belongs  to  the  plane  yz,  the  director  remains  coplanar  at  an  angle 
0  with  respect  to  the  axis  z.  The  angle  0  is  governed  by 

KV20  +  £a  [sin(20)  (|Ey|2  -  |EZ|2  )  +  2Re(EyE*z)  cos(20)  ]  =  0 

with  K  the  Frank  elastic  constant  (for  molecular  deformations).  We  look  for  single-hump 
solitary  solutions  having  a  flat-phase  transverse  profile,  indicating  with  8=8(0)  the  (power 
dependent)  walk-off.  Therefore,  we  fix  the  wavevector  and  allow  the  Poynting  vector  to 
change  its  direction  in  yz,  consistently  with  typical  experimental  set-ups.8,11  Due  to  the  high- 
nonlocality  of  this  medium,19  to  compute  the  soliton  profile  we  can  take  for  8  and  the 
extraordinary  refractive  index  the  value  computed  at  the  beam  peak  (solitons  are  even  owing 
to  symmetry).  In  other  words,  assuming  known  the  distribution  of  0  in  the  self-localization 
regime,  we  can  consider  the  latter  an  inhomogeneous  distribution  of  the  optic  axis,  and  keep 
into  account  that  the  nonlinear  perturbations  are  perceived  as  small  by  the  soliton  because  of 
the  high  nonlocality;  hence,  only  transverse  variations  in  the  refractive  index  are  relevant. 

In  order  to  get  a  simplified  (1+1)D  model  for  nonlinear  optical  propagation,  for  the  x- 
dependence  of  the  optical  perturbation  we  make  the  ansatz  sin  (nx/L)  (the  cell  boundaries  are 
in  x=0,L)  and  consider  beams  propagating  in  the  cell  mid-plane  x=L/2. 

Fig.  3  shows  numerical  solutions  for  0o=tc/4:  due  to  the  reorientational  nonlinearity,  the 
induced  optical  perturbation  saturates  at  high  power,  resulting  in  a  non-monotonic  trend  for 
the  soliton  waist  with  power  [see  Fig.  3(a)]. 
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Fig.  3  Soliton  self-steering  in  NLC.  (a)  Soliton  waist  (dashed  blue  line)  and  walk-off  (solid  green  line)  versus 
power  density  .  (b)  Plots  of  |HX|  and  the  optical  perturbation  A0  in  the  soliton  case,  versus  y  and  the  logarithm  of 
the  power  density  Pd,  with  Pd=62  W/m.  Both  quantities  are  normalized  to  their  peak  value  at  each  power,  (d) 
Beam  evolution  in  yz  computed  by  BPM  when  the  input  profile  corresponds  to  a  solitary  solution  for  Pd  =62 
W/m  (black  lines,  beam  with  the  largest  slope  in  yz)  and  Pd  =1.8  104  W/m  (red  lines),  (c)  and  (d)  show  the 
corresponding  perturbation  in  the  director  distribution  in  the  cell. 

Fig.  3(a)  graphs  also  the  walk-off  8  versus  power,  demonstrating  soliton  self-bending.  Fig.  5 
(c-d)  plot  the  profiles  of  the  field  and  A0  versus  the  transverse  coordinate  y'  and  power 
density  Pd:  the  field  profile  is  very  close  to  Gaussian  due  to  the  high  nonlocality  of  NLC, 
whereas  the  optical  perturbation  A0  is  very  similar  to  the  Green  function  of  the  structure  with 
a  width  slightly  changing  versus  power,  exponentially  decreasing  for  large  |y'|  with  a  slope 
determined  by  L.  Using  a  BPM  propagator  with  a  relaxation  algorithm  to  calculate  the 
director  distribution  in  yz,  we  checked  that  the  calculated  solutions  effectively  correspond  to 
solitons  by  simulating  their  propagation  in  an  NLC  cell.  As  shown  in  Fig.  5  (d),  the  beams 
walk-off  with  invariant  profiles,  confirming  the  validity  of  the  computed  solutions  and  the 
occurrence  of  power-dependent  angular  self-steering. 

4  Conclusions 

The  theoretical  and  numerical  analyses  carried  out  with  reference  to  finite  light  beams  in 
anisotropic  dielectrics  enable  us  to  deal  with  narrow  beams  in  the  non-paraxial  approximation, 
both  in  the  linear  and  in  the  nonlinear  regimes,  detailing  roles  and  contributions  of  the 
uniaxial  response,  the  induced  nonlocality,  the  asymmetry  introduced  by  non-negligible 
longitudinal  field  components  and  the  second  derivatives  in  propagation.  The  detailed 
analysis  of  reorientational  nematic  liquid  crystals  allows  us  to  model  self-steering  of  spatial 
solitons  in  media  with  power-dependent  walk-off  via  reorientation  f  the  optic,  a  regime 
effectively  addressable  in  experimental  samples  with  E7  and  planar  geometries. 
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6.  List  of  figures 

Fig.  1  (a,  d)  Beam  waist  vs  z;  blue  solid  lines  and  red  symbols  refer  to  paraxial  Gaussian 
beams  and  numerical  solutions,  respectively,  (b,  e)  Hx  vs  y/X,  normalized  to  its  peak  and 
computed  for  (from  left  to  right)  zJ  X  =0,  53,  105,  158;  solid  and  dashed  lines  are  HG  profile 
and  numerical  solutions,  respectively,  (c,  f)  Relative  error  between  the  profiles  in  (b,  e).  Top 

and  bottom  graphs  are  computed  for  wo/X  =0.5  and  wo/  X  =2,  respectively.  Here  8//  =2.89, 
s_l=2.25.  Inset:  unit  vectors  s  and  t. 

Fig.  2  (a)  Computed  transverse  profile  of  Re(S)  along  s  and  (b)  relative  error  with  respect  to 
the  paraxial  expression  for  w/X  =0.5;  Re(S)  is  normalized  to  its  peak  for  fixed  z  and  each 
curve  corresponds  to  the  same  z/X  of  Fig.  1.  (c)  Ratio  of  maxima  of  |Et|  and  |ES  |  for  fixed  z 
and  (d)  ratio  between  maxima  of  odd  and  even  parts  of  Re(S)  vs  wfX. 


Fig.  3  Soliton  self-steering  in  NLC.  (a)  Soliton  waist  (dashed  blue  line)  and  walk-off  (solid 
green  line)  versus  power  density,  (b)  Plots  of  |HX|  and  the  optical  perturbation  Dq  in  the 
soliton  case,  versus  y  and  the  logarithm  of  the  power  density  Pd,  with  Pd=62  W/m.  Both 
quantities  are  normalized  to  their  peak  value  at  each  power,  (d)  Beam  evolution  in  yz 
computed  by  BPM  when  the  input  profde  corresponds  to  a  solitary  solution  for  Pd  =62  W/m 
(black  lines,  beam  with  the  largest  slope  in  yz)  and  Pd  =1.8  104  W/m  (red  lines),  (c)  and  (d) 
show  the  corresponding  perturbation  in  the  director  distribution  in  the  cell. 
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SVEA:  slowly  varying  envelope  approximation 
(1+1)D:  one  (transverse)  plus  one  (propagation)  dimensions  in  space 
BPM  :  beam  propagation  method 
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Summary 

In  the  second  period  of  four  months  of  the  Project,  various  aspects  of  light  self- localization  in 
undoped  nematic  liquid  crystals  have  been  examined,  with  particular  emphasis  on  the 
propagation  of  spatial  solitons  and  their  power-dependent  self-steering  via  nonlinear  changes 
in  birefringent  walkoff.  Lor  the  experimental  demonstrations  of  the  results  discussed  hereby 
we  employed  samples  prepared  at  the  University  of  Southampton  (UK)  and  at  the  Institut  Non 
Lineare  de  Nice  (Sophia  Antipolis,  Prance)) 

After  a  brief  introduction,  we  summarize  the  pillars  of  self-localized  beam  propagation  in 
nematic  liquid  crystals,  and  then  report  on  theory,  numerical  predictions  and  experimental 
results  on  transverse  self-acceleration  of  spatial  optical  solitons  in  these  uniaxial  media. 

1  Introduction 

Temporal  and  spatial  self- localization  of  light  are  relevant  topics  in  modem  optics,  for 
their  implications  in  fundamental  physics  as  well  as  their  potential  role  in  new  generations  of 
communication-oriented  all-optical  technologies. 

In  the  spatial  domain,  nonlinear  confinement  of  light  occurs  when  beam  spreading  due  to 
diffraction  is  balanced  by  nonlinear  self-focusing,  giving  rise  to  transversely  finite  beams 
preserving  their  profile  in  propagation,  i.e.  optical  spatial  solitons. [1]  Shape  and  stability  of  a 
soliton  depend  on  the  specific  kind  of  nonlinearity,  and  several  physical  mechanisms  can 
sustain  light  self-trapping,  including  reorientation  [2-3]. 


Dielectrics  with  a  nonlocal  response  are  of  interest  in  light  self-confinement  primarily  for 
the  stabilizing  effects  on  (2+l)D  solitons,  otherwise  unstable  in  pure  Kerr  media.  Nonlocality, 
often  (but  not  exclusively)  quantified  by  the  ratio  between  the  widths  of  the  nonlinear 
perturbation  and  of  the  light-beam  profile,  can  range  from  unity  (local  media)  to  nearly 
infinity  (highly  nonlocal  media).  In  the  latter  limit,  solitons  or  Nematicons  have  been 
observed  in  nematic  liquid  crystals  (NLC)  [2],  Nonlocality  is  also  responsible  for  self¬ 
deflection  of  individual  solitons  via  interactions  with  boundaries,  as  demonstrated  by  us  in 
NLC.  [4] 

We  report  hereby  on  the  investigation,  both  theoretical  and  experimental,  of  nematicon 
self-steering  via  nonlinear  walk-off  and  its  dependence  on  the  initial  orientation  of  the  optic 
axis.  In  Sec.  2  we  model  self-confinement  and  self-steering  of  nematicons  in  the  absence  of 
losses,  solving  the  nonlinear  eigenvalue  problem  for  various  initial  permittivity  tensors  (i.e. 
various  orientations  of  the  optic  axis).  In  Sec.  3  we  report  our  experimental  results  and 
compare  them  with  numerical  simulations  accounting  for  losses. 


2  Methods,  Assumptions,  and  Procedures 

Self-Localized  light  propagation  and  Self-Steering  via  power  dependent  walkoff  in 
nematic  liquid  crystals:  model 

Liquid  crystals  are  mesophase  dielectrics  sharing  properties  of  solids  and  liquids.  The  nematic 
liquid  crystalline  phase  possesses  the  lowest  symmetry,  with  a  high  degree  of  order  in  the 
angular  orientation  of  the  molecules  but  without  positional  order.  We  deal  with  uniaxial 
nematic  liquid  crystals,  i.e.,  with  one  symmetry  axis  or  director  n,  refractive  indices  ne  and  n0 
along  and  normal  to  n,  respectively,  with  ne  <  n0.  We  assume  n  to  lie  in  the  plane  yz:  thus,  the 
angle  0  completely  describes  the  distribution  of  the  director.  It  is  convenient  to  work  with  the 
magnetic  field  H,  with  only  one  non-zero  component  along  n,  rather  than  the  electric  field  E. 
Accounting  for  the  high  nonlocality  of  NLC  stemming  from  intermolecular  links  in  the  liquid 
crystal  state,  nonlinear  optical  propagation  in  the  paraxial  approximation  is  governed  by 

2ikne  dA/dz  +  tan 5  dA/cy  +  Dy  d2 AJdy1  +  d2AJdx2  +  k2An2A  =  0,  (1) 

V  20  +  y|Et|2  sin[2(0  -  8)]  =  0  (2) 

where  k  is  the  vacuum  wavevector,  5  is  the  walkoff  angle,  ne(0)  is  the  extraordinary  refractive 
index,  Dy  is  the  diffraction  coefficient,  A  is  the  slowly  varying  envelope,  An2  accounts  for  the 
light-induced  index  well,  Et  is  the  electric  field  component  orthogonal  to  the  Poynting  vector, 
and  y=£o(n0  -ne  )/(4K),  with  K  the  Frank  elastic  constant  for  molecular  deformations.  Eq.  (2) 
governs  the  director  reorientation  0;  0  affects  the  electromagnetic  propagation  by  modifying 
the  (extraordinary)  refractive  index  Eq.  (1).  At  variance  with  purely  Kerr  media,  the 
relationship  between  refractive  index  changes  and  intensity,  stemming  from  Eq.  (2),  is 
nonlinear.  Specifically,  from  the  sine  term  we  expect  a  saturable  behavior  versus  power, 
strongly  dependent  on  the  initial  conditions.  For  ®  the  optical  perturbation,  we  can  therefore 
set  0=0o+®,  with  0o  a  constant;  determined  by  the  molecular  anchoring  via  chemi-mechanical 
treatments  of  the  interfaces  bounding  the  NLC  layer. 
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Fig.  1.  (a)  Beam  diffraction  in  the  linear  regime,  (b)  Soliton  formation  without  changes  in 
walk-off.  (c)  Soliton  self-steering  with  (high)  input  power,  (d)  Side  view  of  a  typical  planar 
cell. 


We  can  identify  three  different  regimes  for  light  propagation.  At  very  low  powers,  ®  is  small, 
thus  in  the  linear  regime  the  beam  diffracts  and  has  a  Poynting  vector  at  an  angle  5  with 
respect  to  z  [Fig.  1(a)],  When  the  power  increases,  ®  grows  self-confinement  occurs,  but  is 
still  small  with  respect  to  0o  and  does  not  affect  the  walk-off  [Fig.  1(b)]:  this  is  the 
perturbative  nonlinear  regime.  For  further  increases  in  power,  ®  becomes  comparable  with 
0o  and  5  depends  on  power,  i.  e.  nematicon  self-steering  take  place  due  to  nonlinear  variations 
in  walk-off  [Fig.  1(c)]:  this  regime  is  no  longer  perturbative. 

The  system  of  Eqs.  (1-2)  is  a  full  3D  model,  with  no  approximations.  To  simplify  the 
treatment  we  reduce  it  to  an  effective  2D  model  by  considering  a  standard  planar  cell  of 
thickness  L  [Fig.  1(d)]  and  a  nematicon  propagating  in  the  mid-plane  x=L/2,  where  no  forces 
are  exerted  by  the  boundaries.  Using  such  an  effective  2D  model,  looking  for  nematicon 
profiles  versus  0o  and  excitation  power,  we  consider  a  power-dependent  walk-off,  i.e.  solitons 
in  the  form  ©U=u(y)  exp(iknNLz),  ®(x,y,z)=<j>(y,z)  sin(7ix/L),  (|>=  <|>u(y').  where  y’=y-z  tan5.  The 
resulting  model  yields 

2nNLne  k2  u=Dy  <92u/dy'2+k2n2eu  (3) 

(l+tan2S)  52(|)u/5y'2-(7t/L)2(|)u+y(Zo/(ne  cosS)  sin(2((|)u+0o-5))  |u|2=0  (4) 

In  the  perturbative  regime  u,  (|)  and  ne  suffice  to  fully  characterize  a  nematicon,  whereas  in  the 
nonperturbative  regime  8  needs  be  found,  as  well. 

Simplified  solutions  for  the  nematicon  profiles  can  be  found  in  the  highly  nonlocal  limit,  i.e. 
when  a  parabolic  profile  can  be  used  for  the  nonlinear  optical  perturbation. 

The  system  of  Eqs.  (3-4)  was  numerically  solved  with  an  iterative  algorithm  for  nonlinear 

eigenvalue  problems;  at  each  cycle,  the  walkoff  5  was  updated  based  on  the  last  computed  (|)u, 
i.e.  accounting  for  its  nonlinear  variations.  Sample  results  for  the  lowest-order  self-confined 
beam  at  a  wavelength  of  1064nm  are  shown  in  Figs.  (2-3).  The  soliton  |u|2  has  a  Gaussian 
shape,  consistently  with  the  hypothesis  of  high  nonlocality,  whereas  the  nonlinear  index  well 

<|)u  is  proportional  to  exp(-7i|y'/L|).  Nematicon  waist  and  walkoff  versus  0o  and  power  P  are 
plotted  in  Fig.  3  and  agree  well  with  the  theoretical  predictions. 


Fig.  2.  Director  distribution  0o  +  (j)u  (left)  and  corresponding  soliton  profile  \uf  (right)  versus 
y'  for  P=0.1  (solid  red  line)  0.5  (dashed-dotted  blue  line)  and  lOmW  (dashed  black  line). 
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Fig.  3.  From  left  to  right:  plots  of  soliton  waist,  maximum  reorientation  angle  and  walk-off 
versus  P  and  0q. 


3.  Results  and  Discussion 

Soliton  Self-Steering  via  power  dependent  walkoff  in  nematic  liquid  crystals: 
experiments. 

As  pointed  out  in  the  previous  section,  the  model  predicts  soliton  self-steering  owing  to 
nonlinear  changes  in  walk-off.  The  latter  phenomenon  was  recently  addressed  experimentally 
in  Ref.  5  using  a  sample  with  an  orientation  9o=80°  to  optimize  the  observable  beam  deviation 
versus  power:  hereby  we  extend  the  study  to  various  0o 

In  the  experiments  we  used  100pm- thick  planar  cells  filled  with  the  commercial  E7,  with 
birefringence  An=0.2.  Upper  and  lower  glass  interfaces  parallel  to  the  beam  propagation  plane 
yz  were  treated  in  order  to  achieve  various  orientations  of  the  molecular  director  in  the  bulk 
NLC.  We  focused  an  infrared  beam  from  a  Nd:YAG  laser  operating  at  1064nm  to  a  waist  of 
about  5  pm,  coupling  it  into  the  sample  with  wavevector  parallel  to  z  and  in  the  midplane 
x=L/2;  using  the  out-of-plane  scattered  light,  with  a  microscope  objective  lens  and  a  CCD 
camera  we  studied  the  beam  evolution  in  yz  as  a  function  of  input  power  and  initial 
orientation  0o.  When  varying  the  excitation  and  the  orientation,  we  avoided  powers  high 
enough  to  induce  transverse  instabilities  and/or  excessive  beam  fluctuations,  averaging  our 
observations  over  a  large  number  of  acquisitions.  Power  dissipation  in  the  samples  due  to 
scattering  was  the  only  significant  departure  of  our  data  from  the  predictions  illustrated  in  the 
previous  section.  To  correctly  interpret  the  results,  the  beam  power  P  was  expressed  as 
P=Poexp(-2az},  with  Po  the  input  excitation  in  z=0  and  a  the  loss  coefficient,  evaluated  to 


about  5  cm'1  in  our  samples.  Due  to  losses,  in  the  nonperturbative  regime  the  strong 
reorientation  is  gradually  reduced  and,  thus,  the  nonlinear  walkoff  is  non-uniform  along  z, 
giving  rise  to  bent  trajectories  until,  for  large  z  and  back  in  a  perturbative  regime,  they 
become  all  parallel  to  the  linear  Poynting  vector. 

Figs.  (4)  and  (5)  show  the  beam  evolution  and  acquired  trajectories  for  various  0o  and  Po.  For 
0o  =80°,  even  if  the  nonlinearity  was  low  (self-confinement  is  weak)  we  observed  large  self¬ 
steering,  because  a  small  reorientation  produces  a  large  variation  of  walk-off:  at  low 
excitations  the  trajectory  has  no  power  dependence,  while  for  Po  close  to  20mW  the  walk-off 
changes  owing  to  large  reorientation.  For  0o~7O°  the  observed  displacement  is  comparable 
with  the  previous  one,  although  transverse  instabilities  limit  the  observation  to  powers  below 
30mW:  the  trade-off  between  nonlinearity  and  power-dependent  walk-off  results  in  the  largest 
soliton  steering. 

For  0o~45°  the  nonlinearity  increases  towards  its  peak  (the  highest  nonlinear  index  well  for  a 
fixed  power  in  the  perturbative  regime),  but  variations  in  walkoff  decrease:  self-steering 
exhibits  a  relative  minimum  versus  0o,  with  excitation  limited  to  lOmW  by  instabilities. 
Smaller  0o  allow  appreciable  soliton  steering.  From  Fig.  (5)  we  also  infer  that  walk-off 
variations  versus  power  are  positive  for  0o>45°  and  negative  otherwise. 

To  model  the  role  of  losses  on  soliton  propagation  we  used  a  beam  propagator  and  an 
additional  attenuation  term  depending  on  a;  the  nonlinear  changes  in  permittivity  were 
evaluated  from  the  boundary  problem  with  a  relaxation  algorithm;  the  two  solvers  ran 
iteratively  until  convergence  was  reached.  The  beam  trajectories  are  plotted  in  Fig.  (6),  with 
an  optimum  agreement  between  theory  and  experiments  obtained  with  a  fit-coefficient  larger 
than  1  on  soliton  power.  A  comparison  between  measured  and  predicted  power-dependent 
variations  in  walk-off  can  be  carried  out  by  considering  the  self-steering  in  early  stages  of 
propagation  and  at  different  powers  for  all  0o,  namely  Po=2,  5  and  lOmW. 

4  Conclusions 

We  have  investigated  soliton  self-steering  in  undoped  nematic  liquid  crystals  with  various 
initial  orientations  in  the  plane  of  propagation,  i.e.  self-deflection  of  nematicons  owing  to 
power-dependent  reorientation  of  the  optic  axis.  We  introduced  perturbative  and  non¬ 
perturbative  regimes,  the  latter  accounting  for  all-optical  changes  in  walk-off.  We  reduced  the 
complete  3D  system  to  an  equivalent  2D  model  retaining  nonlocality.  In  the  perturbative 
regime,  we  found  two  scalar  parameters  describing  soliton  waist  and  walk-off  versus  0o. 

To  verify  the  predictions,  we  experimentally  investigated  nematicon  propagation  and  self¬ 
steering  in  samples  with  various  0o,  obtaining  an  excellent  agreement  between  data  and 
numerics  even  in  the  presence  of  scattering  losses  which  induce  a  decay  of  the  nonlinear 
response  along  propagation. 
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Fig.  4.  Photographs  of  beam  evolution  for  several  initial  angles  Go  and  power.  Pmax  indicates 
the  maximum  power  before  instabilities  kicked-in.  Pmax  =  [60  30  10  10  10  30]mW  for  0o  = 
[80  70  50  40  30  20]°,  respectively. 


Fig.  5.  Acquired  beam  trajectories  and  their  dependence  on  0o  and  Po.  Panels  (a-f)  correspond 
to  Go  =  80°,  70°,  50°,  40°,  30°  and  20°,  respectively.  Blue,  green  and  red  lines  correspond  to  P 
=  1,5,  lOmW.  The  cyan  line  (where  present)  corresponds  to  the  maximum  usable  excitation, 
Po  =  60m  W  for  0o=  80°  and  Po  =  30m W  for  Go  =  70°  and  20°,  respectively. 
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Fig.  6.  Numerically  computed  soliton  trajectory  in  the  plane  yz  for  various  0o  and  Po— >0mW 
(blue  line),  Po=  lmW  (green  line),  Po=  3mW  (red  line)  and  Po=  5mW  (cyan  line)  versus  0o. 
The  parameters  are  consistent  with  the  experimental  values. 
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List  of  Symbols,  Abbreviations  and  Acronyms 

NL:  nonlinear 

NLC:  nematic  liquid  crystals 

E7:  commercial  mixture  of  nematic  liquid  crystals  (supplied  by  Merck) 
SVEA:  slowly  varying  envelope  approximation 
(1+1)D:  one  (transverse)  plus  one  (propagation)  dimensions  in  space 
BPM:  beam  propagation  method 


Appendix:  confirmation  of  expenditures 

I  confirm  that  expenditures  in  excess  of  US$20,000.00  were  incurred  for  during  the  second 
four  months  (January-April  2011)  of  the  Project  started  last  September  1,  2010,  including 
overhead,  consumables,  equipment,  travel  and  personnel. 

Sincerely, 

Prof.  Gaetano  Assanto 


Grant  #  FA8655-1 0-1 -3010 


NATO-EAORD  Project 

“Re-addressable  Interconnects  with  Light-induced  Waveguides  in  Liquid  Crystals” 

Third  and  Final  Periodic  Report 


Prof.  Gaetano  Assanto 

Nonlinear  Optics  and  OptoElectronics  Lab 

Department  of  Electronic  Engineering,  University  Roma  Tre 
Via  della  Vasca  Navale  84,  00146  Rome-  ITALY 


Table  of  Contents 

Summary 

1 .  Introduction 

2.  Methods,  Assumptions,  and  Procedures 

Self-Localized  beam  propagation,  refractionand  total  internal  reflection  at  a  voltage 
defined  interface  between  two  regions  of  nematic  liquid  crystals 

3.  Results  and  Discussion 

Soliton  waveguide  electro-optic  steering  device:  design  and  numerical  experiments  in 
standard  E7 

4.  Conclusions 

5.  List  of  publications  based  on  the  Project 

6.  References 

7.  List  of  Ligures 

List  of  Symbols,  Abbreviations  and  Acronyms 
Appendix:  confirmation  of  expenditures 
Disclosure  of  inventions 


Summary 

In  the  third  four-month  period  of  the  Project,  various  aspects  of  transverse  acceleration  of 
self-confined  light  beams  in  undoped  nematic  liquid  crystals  have  been  examined,  with 
particular  emphasis  on  the  maximization  of  power-dependent  soliton  self-steering  through  a 
dielectric  interface  defined  by  applied  voltage  in  a  planar  liquid  crystalline  cell.  After  the 
modelling  and  the  detailed  numerical  simulation  of  the  problem,  we  designed  a  specific  set  of 
electrodes  to  maximize  soliton  waveguide  deflection  with  electro-optic  control.  We  expect 
angular  deflections  as  large  as  70  degrees  at  a  single  graded  interface  using  the  standard 
mixture  E7. 

After  a  brief  introduction,  we  describe  the  model,  the  geometry  and  its  advantages;  we  then 
report  on  the  numerical  results  and  expected  performance  on  soliton  steering  upon  refraction 
and  total  internal  reflection. 


1  Introduction 

Owing  to  their  unique  properties  intermediate  between  those  of  liquids  and  of  solids, 
liquid  crystals  (LC)  are  relevant  materials  in  optical  signal  processing;  one  of  their  remarkable 
features  is  their  sensitivity  to  applied  perturbations,  such  as  electric  fields.  The  latter  provides 
large  electro-optic  and  all-optical  responses  and  significant  tunability  through  molecular 
reorientation.  Self-focusing  in  nematic  LC  (NLC),  i.e.  LC  with  a  high  degree  of  orientational 
order  but  no  positional  order  has  allowed  generating  reorientational  spatial  solitons  or 
nematicons,  i.e.,  non-diffracting  self-confined  beams,  at  mW  powers.[l]  Nematicons  are 
light-induced  waveguides  with  nonlocal  features  providing  stability,  robustness  and  long 
range  interactions,  even  with  partially  incoherent  light.  Such  light-induced  channels  can 
confine  and  route  weaker  copolarized  signals  of  arbitrary  wavelengths  and  transmission 
protocols;  they  can  form  a  reconfigurable  fabric  of  interconnections  in  circuits  defined  by 
light  and  with  paths  controlled  by  light  and/or  by  voltage(s). 

We  analyze  a  scheme  for  the  effective  angular  deflection  of  self-guided  beams,  using 
NLC  and  voltage-driven  dielectric  interfaces  to  obtain  the  steering  of  spatial  soliton 
waveguides  over  angles  as  large  as  70°.  We  resort  to  a  double  set  of  comb-patterned 
electrodes  to  achieve  in-plane  nematicon  propagation  and  steering,  by  bias-tuning  the 
effective  anisotropy  and  avoiding  the  detrimental  effects  of  out-of-plane  walk-off. 

Electro-optic  beam  deflectors  with  large  angles  were  previously  demonstrated  in  LC, 
[2,3];  however,  the  use  of  nematicons  with  waist  of  a  few  pm  permits  the  deflection  of 
guided-wave  signals  irrespective  of  their  spectral  placement  and  temporal  format  while 
maximizing  the  number  of  resolvable  output  spots. 


2  Methods,  Assumptions,  and  Procedures 

Self-Localized  beam  propagation,  refractionand  total  internal  reflection  at  a  voltage 
defined  interface  between  two  regions  of  nematic  liquid  crystals 

The  planar  geometry  is  sketched  in  Fig.  1:  a  nematic  liquid  crystal  layer  is  confined 
between  two  glass  slides  parallel  to  plane  yz,  separated  by  L  along  x  and  treated  at  the 
glass/NLC  interfaces  to  ensure  the  uniform  orientation  of  the  optic  axis  n  (i.e.  the  molecular 
director  at  an  angle  0o  with  z.  At  each  interface,  two  interdigitated  sets  of  electrodes  with 
fingers  along  y  and  period  A  permit  the  application  of  low  frequency  electric  fields.  Each 
electrode  pair  on  a  given  interface  has  a  ground  terminal  in  y=0,  such  that  a  voltage  V)  (V2) 
on  the  left  (right)  comb  can  electro-optically  twist  the  NLC  molecular  dipoles  and  therefore 
reorientate  the  optic  axis  in  the  plane  yz  and  throughout  the  layer  thickness  in  the  region 
underneath. 
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Figure  1.  Sketch  of  the  device,  (a)  Front  view;  the  superposed  grey-scale  map  represents 
the  intensity  of  ELf  (black  is  the  highest),  (b)  Top  view  of  the  electrodes,  (c)  Distribution  of 
the  NLC  optic  axis  in  the  cell  mid-plane  x=L/2  (bottom)  and  corresponding  0  versus  y  (top). 
Inset:  definition  of  0.  The  red  arrows  indicate  the  orientation  of  the  molecular  director. 


For  A  «  L,  the  dominant  component  of  the  applied  electric  field  is  along  y,  thus  n  tends 
to  rotate  in  the  plane  yz  and  align  with  the  field  polarization.  To  first  approximation,  at  each 
NLC/glass  interface  ELf(x=0)=Elf(x=L)=0  under  the  electrodes  and  ELf(x=0)= 
Elf(x=L)=2Vi.2/  A  in  between  them;  the  Elf  distribution  for  V2>Vi  is  represented  in  Fig.  1(a). 
Owing  to  the  NLC  electro-optic  response  to  Elf  the  optic  axis  n  will  be  twisted  everywhere  in 
the  same  way.  Moreover,  the  nonlocal  character  of  the  elastic  response,  in  the  limit  A  «  L, 
will  smooth  out  the  spatial  inhomogeneities  of  Elf,  making  the  bias-driven  reorientation 
voltage-tuned  in  the  two  regions  across  y=0  and  uniform  through  the  cell  thickness. 
Therefore,  applied  biases  Vi,2  give  rise  to  a  graded  interface  along  y  [see  Fig.  1(c)],  of  width  d 
depending  on  L  and  with  an  electrically  controlled  change  in  director  angle. 


Figure  2.  (a)  Reorientation  angle  in  the  cell  mid-plane  x=L/2  for  Vi=V2  versus  applied 
bias:  each  label  indicates  the  corresponding  background  angle  0o;  (b)  color  map  of  0  in  the 
plane  xy  and  (c)  maximum  change  in  0  (x=L/2,y)  versus  applied  biases  Vi  and  V2.  (d-f):  0 
distribution  versus  y  in  x=L/2  for  various  V2.  The  labels  on  each  line  indicate  the 
corresponding  Vi  The  bias  values  refer  to  a  thickness  L=  100pm. 

We  numerically  solved  the  NLC  reorientational  equation 

926lf  920lf  .  \  1 7-1  1 2  n 

— I — — \~1lf  sm  (2 9lf)  I  Elf  \ 1  =  0> 

with  0lf  the  electro-optic  reorientation,  Ylf=Aslf/  [2K]  with  AsLf  the  low-frequency 
dielectric  anisotropy  and  K  the  elastic  constant.  We  chose  the  parameters  of  the  commercial 
NLC  mixture  E7,  with  Ylf=5.4V~2;  we  also  employed  g=3.3.  To  trade  off  voltage-sensitivity 
(the  Coulombian  torque  is  maximum  for  0o=45°)  and  the  0  variation  across  the  interface,  we 
set  0o=lO°  [see  reorientation  curves  in  Fig.  2a].  As  expected,  0  is  constant  near  the  mid-plane 
x=L/2  and  undergoes  a  negative  (positive)  variation  A0  along  y  for  Vi>V2  (Vi<V2)  [Fig.  2c], 
with  a  quasi-linear  trend  over  a  transition  length  d~L,  independent  on  bias. 

Having  established  the  link  between  voltages  and  0  distribution,  i.e.  refractive  index  for 
extraordinarily  polarized  light,  we  turn  to  nonlinear  propagation  of  beams  launched  in  y<0 
(where  the  index  gradient  is  negligible)  with  wave  vector  k  at  an  angle  a  with  z  in  the  plane 
x=L/2  and  extraordinary  field  polarized  along  y.  Such  excitation  can  reorientate  the  optic  axis 
in  the  plane  yz,  i.e.  ensuring  in-plane  director  rotations  and  corresponding  in-plane  changes  of 
walk-off.  We  stress  that  the  spatial  soliton  trajectory  after  the  interface  does  not  depend  on  the 
profile  of  the  transition  region  but  simply  on  A0  ,  the  specific  distribution  of  0Lf  only 
affecting  the  lateral  displacement  of  the  beam. 

Nonlinear  reorientation  in  yz  and  self-confined  nematicon  dynamics  can  be  modelled  in 
2D  by 
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with  A  the  magnetic  field  along  x, 

Y  =  [e0/(4if)]  («|  -  )  [Zo  /(no  cos  <5)]2  (Z0  is  the  vacuum  impedance), 

y/  the  (nonlinear)  perturbation  on  0  due  to  light  (i.e.,  0  =0lf  +  Y), 

Dy  the  diffraction  coefficient,  5  the  walk-off,  ko  the  vacuum  wavenumber,  no  the 
extraordinary  index  and 

Ani  —  n%  {9)  —  n%  the  overall  refractive  index  change. 

The  nematicon  trajectory  can  be  calculated  using  the  Ehrenfest's  theorem. 


3.  Results  and  Discussion 

Soliton  waveguide  electro-optic  steering  device:  design  and  numerical 
experiments  in  standard  E7 

We  assumed  a  linear  profile  for  0Lf  across  d,  from  y=0  to  y=L=100  pm.  Fig.  3  a  shows 
the  evolution  of  a  Gaussian  beam  launched  in  y=-50  pm  $  with  kllz.  For  0i>  02  the  refractive 
gradient  is  negative,  thus  the  soliton  wave  vector  bends  towards  y<$;  to  achieve  TIR,  the 
change  in  k  has  to  ensure  a  Poynting  vector  S  leaving  the  interface  with  a  negative  y 
component.  Noteworthy,  reflection  and  incidence  angles  are  generally  different  owing  to 
anisotropy.  For  0i<  02,  conversely,  the  beam  can  always  overcome  the  index  barrier  and  be 
refracted. 

Fig.  3b  graphs  the  calculated  soliton  trajectories  versus  input  power:  the  visible  power 
dependence  is  due  to  all-optical  changes  in  walk-off.  For  excitation  with  kHz  the  angle  (3  of 
the  Poynting  vector  with  z  away  from  the  interface,  the  overall  steering  (from  refraction  to 
TIR)  amounts  to  about  40°  [Fig.  3c-d],  with  excellent  agreement  between  the  results  from  a 
beam  propagator  (BPM)  and  the  predictions  of  the  geometric  optics  approximation  (GOA)  for 
P  up  to  20°.  Even  larger  deflections  are  obtained  by  launching  solitons  with  phase  fronts 
slanted  with  respect  to  z,  i.e.  nonzero  a. 
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Figure  3.  (a)  Nematicon  evolution  for  P=3mW;  dashed  and  dotted  white  lines  indicate  the 
interface  borders  and  the  beam  trajectory  when  0Lf  =10°,  respectively,  (b)  Left:  trajectories 
upon  TIR  for  P=0.5mW  (no  symbols)  and  P=3mW  (squares);  right:  dashed  (solid)  lines 
correspond  to  02=90°  (10°).  Soliton  output  slope  (3  versus  0i  for  02=10°  (c)  and  versus  02  for 
01=10°  (d).  The  wavelength  is  1064nm,  the  input  waist  5pm  and  L=100  pm. 


Figure  4.  (a)  TIR  region  vs  0i  and  02.  Angle  (3  versus  02  for  a=0  (b),  25°  (c)  and  70°  (d); 
the  labels  indicate  the  corresponding  0i.  (e)  Maximum  angular  deflection  of  wave  vector  and 
Poynting  vector  versus  a. 

For  a  fixed  a,  the  crucial  requisite  to  achieve  large  steering  is  the  occurrence  of  TIR. 
Defining  a™  the  maximum  angle  of  incidence  resulting  in  TIR,  a™  is  plotted  in  Fig.  4a  \  for 
given  0i  2,  with  a  maximum  oitir  approximately  equal  to  25.1°;  above  such  curves  TIR  does 
not  take  place.  The  output  angle  (3  versus  02  and  0i  is  graphed  in  Fig.  4b-d:  the  largest 
deflections  are  obtained  for  a<max(aiiR)  [for  a  fixed  0i  the  solitons  evolve  from  TIR  to 


refraction  as  02  varies,  see  Fig.  4(b-c]  rather  than  for  a>max(axiR)  [where  only  refraction 
occurs,  see  Fig.  4d],  with  a  maximum  steering  of  60°  for  a=max(axiR)  and  0i=9O°.  By 
varying  0i  as  well,  an  overall  steering  angle  of  70°  is  achieved  for  max(axiR)  [Fig.  4e],  taking 
advantage  of  the  dependence  of  the  refraction  angle  from  walk-off. 

Such  voltage-controlled  angular  deviations,  to  be  achieved  with  biases  of  a  few  Volt,  are 
among  the  largest  observed  at  a  single  dielectric  interface  and  well  above  the  record  value  of 
$40°  reported  for  spatial  solitons.[4]  Larger  steering  could  be  engineered  by  adopting  more 
complex  or  multiple  interfaces. 


4  Conclusions 

We  introduce  a  novel  concept  for  optical  signal  steering  using  spatial  solitons  and  a  voltage- 
tunable  interface  in  nematic  liquid  crystals.  The  design  ensures  in-plane  self-confined  and 
signal  guiding  propagation  and  deflection  over  angles  as  large  as  70°  by  employing  modest 
voltages,  enabling  the  realization  of  readdressable  interconnections  and  spatial  demultiplexers 
with  a  wide  angular  span  while  guaranteeing  the  transparent  transmission  of  one  or  several 
signals  with  arbitrary  wavelengths  and  modulations. 


5  List  of  Publications  based  on  the  Project  in  International  (ISI)  Journals 

-  A.  Alberucci  and  G.  Assanto:  “Nematicons  beyond  the  perturbative  regime“,  Opt.  Lett.  35 
(15),  2520-2522  (2010) 

-  A.  Piccardi,  A.  Alberucci  and  G.  Assanto:  “Power-dependent  nematicon  steering  via  walk- 
off’,  J.  Opt.  Soc.  Am.  B  27  (1 1),  2398-2404  (2010) 

-  A.  Alberucci  and  G.  Assanto:  “Non  Paraxial  Solitary  waves  in  anisotropic  dielectrics”,  Phys. 
Rev.  A  83,  033822  (2011) 

-  A.  Alberucci  and  G.  Assanto:  “On  Beam  Propagation  in  Anisotropic  Media:  one¬ 
dimensional  analysis”,  Opt.  Lett.  36  (3),  334-336  (2011) 

-  A.  Alberucci  and  G.  Assanto:  “Non-paraxial  (1+1)D  spatial  solitons  in  uniaxial  media”,  Opt. 
Lett.  36(2),  193-195  (2011) 

-  A.  Piccardi,  A.  Alberucci,  N.  Tabiryan  and  G.  Assanto:  “Dark  Nematicons”,  Opt.  Lett.  36 
(8),  1456-1458  (2011) 

-  W.-P.  Zhong,  M.  Belie,  G.  Assanto,  B.  A.  Malomed,  and  T.  Huang:  “Self-trapping  of  scalar 
and  vector  dipole  solitons  in  Kerr  media”,  Phys.  Rev.  A  83,  043833  (2011) 

-  C.  P.  Jisha,  A.  Alberucci,  R.-K.  Lee,  and  G.  Assanto:  “Wave-Particle  Duality  of  Optical 
Solitons  interacting  with  defects”,  Opt.  Lett.  36  (10),  1848-1850  (2011) 

-  M.  Kwasny,  A.  Piccardi,  A.  Alberucci,  M.  Peccianti,  M.  Kaczmarek,  M.  Karpierz  and  G. 
Assanto  “Nematicon-nematicon  interactions  in  a  medium  with  tunable  nonlinearity  and  fixed 
nonlocality”,  Opt.  Lett.  36  (13),  2566-2568  (201 1) 

-  R.  Barboza,  A.  Alberucci  and  G.  Assanto:  “Large  electro-optic  beam  steering  with 
Nematicons”,  Opt.  Lett.  36  (14),  2611-2613  (2011) 

-  G.  Assanto,  A.  Alberucci  and  A.  Piccardi:  “Transverse  self-acceleration  of  Nematicons:  can 
a  self-confined  beam  change  its  own  path?”,  J.  Nonl.  Opt.  Phys.  Mat.,  submitted  (in  press) 


6  References 

1  G.  Assanto  and  M.  Karpierz,  Liq.  Cryst.  36,  1161  (2009) 

2  S.  Khan  and  N.  Riza,  Opt.  Express  12,  868  (2004) 

3  S.  Davis,  G.  Farca,  S.  Rommel,  A.  Martin,  and  M.  Anderson,  Proc.  SPIE  6971,  69710G 
(2008). 

4  M.  Peccianti,  A.  Dyadyusha,  M.  Kaczmarek,  and  G.  Assanto,  Nat.  Phys.  2,  737  (2006). 

7.  List  of  figures 

Figure  1.  Sketch  of  the  device,  (a)  Front  view;  the  superposed  grey-scale  map  represents  the 
intensity  of  Elf  (black  is  the  highest),  (b)  Top  view  of  the  electrodes,  (c)  Distribution  of  the 
NLC  optic  axis  in  the  cell  mid-plane  x=L/2  (bottom)  and  corresponding  0  versus  y  (top). 
Inset:  definition  of  0.  The  red  arrows  indicate  the  orientation  of  the  molecular  director. 

Figure  2.  (a)  Reorientation  angle  in  the  cell  mid-plane  x=L/2  for  Vi=V2  versus  applied  bias: 
each  label  indicates  the  corresponding  background  angle  Go;  (b)  color  map  of  0  in  the  plane  xy 
and  (c)  maximum  change  in  0  (x=L/2,y)  versus  applied  biases  V)  and  V2.  (d-f):  0  distribution 
versus  y  in  x=L/2  for  various  V2.  The  labels  on  each  line  indicate  the  corresponding  V)  The 
bias  values  refer  to  a  thickness  L=  100pm. 

Figure  3.  (a)  Nematicon  evolution  for  P=3mW;  dashed  and  dotted  white  lines  indicate  the 
interface  borders  and  the  beam  trajectory  when  0Lf  =10°,  respectively,  (b)  Left:  trajectories 
upon  TIR  for  P=0.5mW  (no  symbols)  and  P=3mW  (squares);  right:  dashed  (solid)  lines 
correspond  to  02=90°  (10°).  Soliton  output  slope  p  versus  0i  for  02=10°  (c)  and  versus  02  for 
01=10°  (d).  The  wavelength  is  1064nm,  the  input  waist  5pm  and  L=100  pm. 

Figure  4.  (a)  TIR  region  vs  0i  and  02.  Angle  p  versus  02  for  a=0  (b),  25°  (c)  and  70°  (d);  the 
labels  indicate  the  corresponding  0i.  (e)  Maximum  angular  deflection  of  wave  vector  and 
Poynting  vector  versus  a. 


List  of  Symbols,  Abbreviations  and  Acronyms 

NL:  nonlinear 

LC:  liquid  crystals 

NLC:  nematic  liquid  crystals 

E7:  commercial  mixture  of  nematic  liquid  crystals  (supplied  by  Merck) 

BPM:  beam  propagation  method 
GOA:  Geometric  optics  approximation 
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